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Abstract 

We begin to study in tliis paper orbital and asymptotic stability of standing waves for a model of Schrodinger 
■ equation with concentrated nonlinearity in dimension three. The nonlinearity is obtained considering a point 
(or contact) interaction with strength a, which consists of a singular perturbation of the laplacian described 
by a selfadjoint operator Ha, where the strength a depends on the wavefunction: iii = H^u, a — a{u). If q 
is the so-called charge of the domain element u, i.e. the coefficient of its singular part, we let the strength a 
depend on u according to the law a — -~i'\q\'' , with i/ > 0. This characterizes the model as a focusing NLS with 
concentrated nonlinearity of power type. For such a model we prove the existence of standing waves of the form 
u{t) — e*"*<i>t^, which are orbitally stable in the range a e (0, 1), and orbitally unstable for cr > 1. Moreover, we 
show that for a £ (Oj every standing wave is asymptotically stable in the following sense. Choosing initial 
data close to the stationary state in the energy norm, and belonging to a natural weighted space which allows 
dispersive estimates, the following resolution holds: u{t) — 6*"=°*$!^^ + Ut * ipoo + r^o, as t — > +oo, where 
U is the free Schrodinger propagator, Uoo > and V'oo, ''oo G L^(R'^) with ||roo||L2 = 0(t^'''/'*) as i ^ +oo. 
Notice that in the present model the admitted nonlinearity for which asymptotic stability of solitons is proved 
is subcritical. 
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I. INTRODUCTION 



In this paper we begin a systematic analysis of the stability of solitary waves for a nonlinear 
Schrodinger equation with a nonlinearity concentrated in space dimension three. In particular, we 
show that the standing waves of the model are asymptotically stable in the sense that the evolution of 
the system in a neighbourhood of a standing solution admits a soliton resolution expansion: at large 
times, the evolution decomposes as the sum of a standing wave (possibly with different parameters 
from those of the reference initial soliton), a free linear wave, and a small remainder with a spatial 
decay stronger than the linear dispersive one. 

An analogous study concerning the NLS equation with a concentrated nonlinearity in dimension 
one was given in [4| and 



22( 1 . These papers have been a source of inspiration for the present work, 
in particular for what concerns the general scheme of analysis and for some proofs. However, the 
one and the three-dimensional models are different, in particular the latter is strongly singular and 
its energy space is not contained in i7^(IR^). This fact prevents us from following step by step the 
techniques and the results of the cited papers; in particular, no formal manipulations with delta 
distributions are possible, and the full definition of a delta interaction as a point perturbation of the 
laplacian is needed in the analysis. We shall comment on that along the paper. 
We start by giving a presentation of the model. According to pj], we construct a Schrodinger equation 
with concentrated nonlinearities in dimension three by starting from the standard three-dimensional 

n 

linear Schrodinger operator with a so-called point or delta interaction ([3]). Point interactions are 
widely used in Quantum Mechanics as models of contact or zero-range interactions and they are 
intended to describe strongly concentrated potentials at a point. In order to rigorously define a delta 
interaction located at the origin of M'^ we first consider the Laplacian restricted to the set Cg"(M'^\{0}) 
and obtain a symmetric non selfadjoint operator with deficiency indices (1, 1). Second, by the classical 
Von Neumann-Krejn theory there exists a one-parameter family of selfadjoint extensions, which we 
denote by Ha- The operator Ha is defined on the domain 

D{Ha) = {uG L^{R^) : u{x) = </)(x) + qGo{x) with G L^M"^) , V0 G L^{M.^) , A(j) G L^{R^), 

q £ C, lim(u(x) — qGo{x)) = aq}, 
where Go is the Green's function of the Laplacian in three dimensions, i.e. 
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and the action is given by H^uix) = — A</>(x), x G M^. To summarize, any element of the domain 
decomposes in a regular part (j) and a singular (Coulombian) part; the coefficient q of the singular 
part is conventionally called charge^ and the boundary condition imposes a relation between the 
charge and the value of the regular part at the origin depending on the so-called strength a of the 
point interaction, which is the parameter that fixes the selfadjoint extension. 

An alternative equivalent and perhaps more direct construction, which better justifies the interpre- 
tation and the physical meaning of H^^ can be given by defining as a suitable scaling limit (in 
norm resolvent sense) of a family of Schrodinger operators of the form — A -|- VJ,, where is a short 
range potential that approximates a delta distribution as e — )• 0. Performing such limit requires a 
rescaling procedure in order to yield a non-trivial result, and the parameter a appearing in the above 
definition characterizes the particular selfadjoint extension and is related to zero energy resonances 
of the approximating operators. For details and futher information see jsl. 

Whatever the definition given to the operator is, we recall that, for q > (repulsive delta 
interaction). Ha is positive and its spectrum is purely absolutely continuous and coincides with 
[0,+oo), while for a < (attractive delta interaction) an isolated simple negative eigenvalue A = 
— (47ra)'^ appears, corresponding to a bound state. A second property relevant to the physical 
interpretation of the model and related to the value of a is that the scattering length of a delta 
interaction of strength a is given by —{A-na)^^ . The closed and lower bounded quadratic form 
associated to Hn is 



y = G L2(m3) : u{x) = (t)u{x) + quGoix), with G ^^(K^), V(^„ G L\R^), G C}, (3) 



Note that for a generic element u of the form domain the charge q and its regular part (j) are 
independent of each other. Note also that the energy domain is strictly larger than H^{W^). So, the 
linear problem cannot be considered as a small perturbation of the standard free problem in the sense 
of the quadratic forms (at variance with the one-dimensional case). An equivalent representation of 




(2) 



defined on the domain of finite energy states 



which is a Hilbert space endowed with the norm 



u\\v = l|V</'n||L2 + \qu 



(4) 
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the energy space is obtained, fixed A > 0, by 

V=i^u = cl)x+ qGx, withcbx G H'{R^), q€C, Ga(x) = , (5) 

and one can define an equivalent energy norm by 

\W\\v = \NMh + \qf, yueV. 

Notice that G\ £ L^(R'^) and (px £ H^{W^), while in the representation ^ of the regular part was 
just in the homogeneous Sobolev space D^(]R^). 

Following [ij, the nonlinear model can be defined by allowing the strength a to depend on u as 
a{u) = —u\q\^'^^ with > 0, cr > 0, so that 

D{Ha{u)) = {ue L2(m3) : u{x) = + qGo{x) with (/) G Hf^^{M.^), A0 G ^^(M^), 

g G C, lim(u(x) — qCoix)) = —I'lql'^'^q}, 

and H^i^y^u = —Ac/). In the following sections, we often omit the notation -ffQ,(u) in favour of Ha if no 
risk of confusion exists between the linear and the nonlinear operator. We stress that the nonlinarity 
we are considering is focusing. It can be interpreted as modeling the action of a defect in a medium 
which exerts a nonlinear response to the propagation. We remark that a more general definition of 
concentrated nonlinearities (with applications to the case of the wave equation) is given in [ill]. 
We consider the evolution generated by the nonlinear operator Ha[u)^ i-S- 

du 

I— = Ha{u)U. (6) 

In the present literature, there is some physical and numerical analysis of Schrodinger dynamics in 
the presence of nonlinear defects, mainly focused on the milder one-dimensional case ( [2^ , j^] j |3] ) ■ 
The more technical construction of the three-dimensional problem has hindered extended modelistic 
study, numerical work as well as rigorous analysis. Moreover, a certain amount of literature is devoted 
to NLSE with nonhomogeneous (i.e. x-dependent and decaying) nonlinearities, yet with a relatively 
low decay at infinity (see Q, [3] and references therein) . 

Local (for any o" > 0) and global (for o" < 1) well-posedness of the Cauchy problem associated to the 
nonlinear Schrodinger equation ([6]) in the space V have been established in S and 0]. In particular, 
^ admits two conserved quantities called mass and energy, defined as 

M{u{t)) = \\u{t)\\l,, E{u{t)) = - s^kWI'^^^'- 



In Section [IT] we prove that equation ^ admits standing waves, i.e. solutions of the form u{x,t) = 
e*'^*<I>t^(x) , where the profile or amplitude up to a phase factor e*^ is given by 

yAiTu J A7r\x\ 

The set of standing waves is called the solitary manifold A4, and the main concern of this paper 
consists in the study of the large-time evolution of initial data in the vicinity of 7W. A first result 
concerns stability and instability of standing waves. Stability has to be intended as orbital stability, 
i.e. Lyapunov stability up to symmetries of the equation, in this case up to gauge (U (1)) invariance. 
The orbit of is then ^{^co) = {e'^^Luix), 9 G M}. Thus, by definition, the state <I>(^ is orbitally 
stable if for every e > there exists 5 > such that 

d(V'(0),^($^)) < 5 ^ diijit),ffi<^>^)) <e Vt>0 

where d{ip, ^{^ui)) = infug^($^) ||^/) — tiHy. A stationary state is said to be unstable if it is not stable. 

Then, we have the following result, proved in Section Hill 

Theorem (Orbital Stability) Let us consider ([6]). Then, for every cj > 0, 

(a) if < (7 < 1, then the state is orbitally stable 

(b) if o" > 1, then is orbitally unstable. 

The result directly follows from Weinstein 30] and Grillakis-Shatah-Strauss theory for the 
case o" 7^ 1, while for the case a = 1 the pseudoconformal invariance of the equation gives the 
instability by blow-up. 

The core of the paper is devoted to the study of the asymptotic stability of the family of stationary 
states. Asymptotic stability means, loosely speaking, that the solution u{t) corresponding to an 
initial datum u{0) close to the family of orbits, approaches some element of the family of orbits as 
t — )• 00. The analysis makes use of the representation 

^(t,x) = e^eW(<l>^(i)(x)+x(t,x)), (8) 

where Q{t) = fQUj{s)ds + 7(t), and 7(t) is a suitable phase. The solution is represented at every 
time as a modulated solitary wave, with time dependent parameters, up to a fluctuating remainder 
X which has to be controlled. Asymptotic stability of the family of standing waves means that the 
modulating parameters u}{t) and j{t) have a limit as t — t- 00, and the fluctuation x is in some sense 
a small and decaying dispersive correction; the radiation damping through dispersion is responsible 
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for the "dissipative" asymptotic behaviour of the solution u around the family of relative equilibria 
ff{^u}) ■ Notice that, however, in general the solution does not converge to the solitary wave to which 
it was close initially. 



The subject of asymptotic stability of solitary waves was pioneered by Soffer and Weinstein ([25l|. 
26|), and Buslaev and Perelman ([a], ^), who developed the main strategies and techniques, nowa- 
days classical; a more recent presentation is contained in ^2]- Many relevant later contributions 
refining and enlarging the hypotheses in the original papers, as well as concerning the kind of initial 
admitted data and nonlinearities, are contained in Q, 0, 21, 28, 29]. According to this con- 



solidated analysis, one must preliminarily indagate the spectrum of the linearization of equation ([6]) 
around the solitary solution. Writing u = e*'^*(<l>^ + R) and identifying R with the vector of its real 
and imaginary part, we obtain that it satisfies the canonical system 



jdR 
'dt 



Ha^+OJ 

Hc,.^+u 



R = DR 



where H^^ are (linear) delta interaction hamiltonian operators with fixed strength Uj that depend 
on the stationary state <1>^ (through its charge) and on the parameters of the model v, a (see eq. 
(jl7p ). So the dynamics of the linearization of the NLSE around the standing wave is controlled 
by the nonselfadjoint (Hamiltonian) matrix operator L = JD. The explicit characterization of the 
spectrum of the linearization L is possible due to the detailed knowledge of the properties of operators 
Ha^. Such feature is unfrequent and allows to avoid further spectral assumptions. The complete 
result is given in Section IIV^ Theorem IIV.2I Here it is sufficient to recall that in this paper we study 
asymptotic stability of standing waves in the range a G (0, l/\/2) only, which corresponds to L having 
no eigenvalues different from zero and no resonances at the threshold of the essential spectrum. A 
forthcoming paper will treat the case a G [l/\/2, 1), where two simple eigenvalues ±i2a\/l — a'^u 
appear. 

Let us notice that the representation ([8]) amounts in fact to a change of coordinates from the original 
global u to the new set {uj, 7, x}, with a finite dimensional component given by {oj, 7}, that describes 
the solitary manifold and an infinite dimensional one described by x- However, the representation is 
not unique, because any choice of a;,7 gives a corresponding choice of x such that u given by ([8]) is 
a solution of ([6]); so one has to restrict in some way the behaviour of the new parameters {a;, 7,%} 
of the solution. To this end, we exploit the fact that the solitary manifold can be naturally endowed 
with a symplectic structure (see Section III Ap and it turns out that its tangent space coincides 



with the generahzed kernel of the hnearization L. The generahzed kernel is in turn non trivial, so the 
propagator e~*^ has a component growing in time. A parametrization of the running approximate 
solitary wave in the neighborhood of the solitary manifold suitable for asymptotic analysis is hence 
obtained through a symplectic splitting in a component along the solitary manifold and a component 
transversal (symplectically orthogonal) to it. Requiring that the infinite dimensional component x 
is purely transversal, i.e. projects to zero on the directions of the discrete spectrum, here reduced 
to the generalized kernel of the linearization, provides the set of the so called modulation (coupled) 
equations for the parameters oj{t) and 7(t), as well as a corresponding partial differential equation 
for X (see [isj for an enlightening description of the symplectic projection method). The goal is to 
establish the asymptotic behaviour of the solutions to the modulation equations with a simultaneous 
control of the decay of the nonlinear part x, through the so-called majorant's method (see ^7]). 
The main result of this paper is the following, and it is proven in Section IVIIi 
Theorem (Asymptotic stability) 

Assume a G (0, l/\/2). Let u{t) G C(M+, V) be a solution to equation ^ with n(0) = G ^nL^ 
and denote d = \\uq — e^^^'^I'tjo llvnLi, i for some ujq > {) and 6q G M. Then, if d is sufficiently small, the 
solution u{t) can be decomposed as follows 

u{t) = e''^-*^^^ + * Voo + r-oo , (9) 

where ojoo > and ipoo, foo G i>^(K^), with ||r-oo||L2 = 0(t^^/^) as t — )• +oo. 

In the previous theorem, is defined in Section HV Bl and is a weighted space of integrable functions. 
The weight guarantees the validity of the dispersive estimates needed in order to control the decay 
of the transversal evolution, and it seems at present unavoidable in view of the singularity of finite 
energy states. Moreover, it imposes a certain localization on the the admitted initial data, which 
seems to be a technical requirement. 

Concerning the treatment of the modulation equations, one of the main additional difficulties with 
respect to standard models, and in particular with the case of concentrated nonlinearities in one 
dimension treated in Im and 



22] , is that the equations controlling the evolution of the transversal part 
X have domains that change with time. This fact forced us to make use of the variational formulation 
(i.e. in terms of quadratic forms) instead of the traditional strong formulation (i.e. in terms of 
operators and their domains) . The same problem propagates to the proof of the asymptotics given 
in the above theorem. A last remark concerns the seemingly anomalous value of the nonlinearities 
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where asymptotic stability is proven; this because in the typical situations, when standard NLS with 
or without potential is treatedjjt is difficult to have information about subcritical nonlinearities (but 
see the notably exception in [l3|), and in particular pure power. On the other hand, the present 
model corresponds to an inhomogeneous (space dependent and strongly singular) nonlinearity; this 
seems to indicate that the analysis of specific models can give results not accessible to general theory, 
at least at present. 

Acknowledgments The authors are grateful to Gianfausto Dell'Antonio and Galina Perelman for 
several discussions. R. A. is partially supported by the PRIN2009 grant Critical Point Theory and 
Perturbative Methods for Nonlinear Differential Equations. 
This paper is part of the Ph.D. thesis of C. O. 



II. PRELIMINARIES 



A. Hamiltonian structure 



We consider ^^(M^C) real Hilbert space endowed with the scalar product 

{u,v)j^2=^ / uvdx= / (Kf 5fiu + 9u9u)dx. (10) 



It is sometimes convenient to shift from the complex valued representation of u to the vector real 
valued one through the identification u = ^fiu + iQu i— )• (5Rm, 9n) = (ni,n2). As a consequence, 
H^{M.^,C) = H^(W^,M.'^), while multiplication by i is equivalent to multipication by the matrix — J, 
where 



J 



1 
-1 



(11) 



The space L (M'^) is also a symplectic manifold when endowed with the symplectic form 

Q{u,v)=Q / uv dx = {^vQu — '^v^u)dx = / {u2Vi — uiV2)dx. (12) 



Along the paper we often shift between real and complex representation when no ambiguity occurs. 
In our model the Hamiltonian functional coincides with the total energy and it is given by 

E{uit)) = ^\\Vm,-^^\q\^'^+\ u = cP + qGoeV. (13) 



Correspondingly, the NLSE ^ takes the hamiltonian form 

du 



J E {u{t)) . (14) 



B. Standing waves 



Standing waves are solutions of the form u{x,t) = e*'^*<I>tj(x). It immediately follows that if a 
standing wave exists, then the amplitude satisfies the nonlinear equation 

H^^^)<^>^ + u;<^^ = 0. (15) 

Proposition II. 1. Standing waves for equation ([6]) exist if and only if > 0. In such a case the 
set of solitary waves is given by the two-dimensional manifold 

M = { e^® , w > , G G [0,27r)} , (16) 

where the function 



^Airu J Att\x\ 

and the parameters co and play the role of local coordinates. 

Proof. Recall that the function Gq defined in ([T]) satisfies the equation — AGq = S 

where 6 is the Dirac's delta distribution centred at x = 0. Hence, for x ^ equation (jl5p is equivalent 

to — A$^(x) +Lo^uj{x) = 0. Consider the corresponding equation in spherical coordinates, namely 

d'^u 2 du 1 d'^u cos (j) du 1 d'^u ^ 

J. Qj. J.2 Qfjj2 J.2 gjj^ ^2 gj]^2 ^ QQ2 ' 

and exploit the spherical harmonics expansion of the solution u{r, 6, (p) = J2tS} 'll\=~i *^)' 
where Y^j denotes the set of spherical harmonics which is an orthonormal basis of L^([0,7r] x 
[0,27r], sin 6l(i6'#). Since 

1 d'^Yi, COS0 dYi^ 1 d'^Yi 



^ . 2 1 aa-? = -^^IJ-: foi' some A G C, 



J.2 Q^2 J.2 gjj^ Q^p jz. gjj^. 

one has that A belongs to the set {A; := lif + 1), / G N}, and so the functions u^j solves —u'l-{r) 
^u'i j{r) + (uj- uij{r) = 0. Then, from formula 8.491.6 in Q, 

%Kr) = ^Zr^^iV^r), 



where Zy is a generic Bessel's function. By the asymptotic expansions 8.443 and 8.451.1 in [iSJ] one 
immediately has that if A 7^ 0, then Uj^i cannot belong to L^(]R+, r^dr). Hence, we fix A = and 
denote $a;(x) = r = Thus u has to be a square-integrable solution of u"(r) — u}u{r) = 0, 
and 

Att\x\ 

for some q € C and u > and c = 0. Furthermore, by the boundary condition in the definition 
of D{Ha(u)), i-e. liioa^^Q{^^{x) - qGo{x)) = a{^^)q, one gets = -i^\q\'^''q, and supposing 

1/ / one obtains \q\'^'^ = This requires > 0, so 

\47ri/ J 47r|x| 

which, up to a phase factor, gives the stated result. In the case = 0, from boundary condition we 
get (7 = or a; = 0. If q = , then the function u vanishes. If a; = 0, then one has u{x) = , which 
is the resonance function of the delta interaction with vanishing strength, but it is not an element 
of the operator domain, and it does not solve the stationary equation p5|) . So for = standing 
waves do not exist. □ 



In the following we denote q^^ = ^3^^ 



1 

2a 



C. Linearization of around 

The linearization of equation ^ around a stationary solution is not completely obvious, due to 
the fact that the nonlinearity is embodied in the domain of the operator Ha(^u) C'l^d not in the action 
of the operator itself. Nevertheless, we can consider the Hamiltonian associated to equation ([6]) given 
by formula ()14p and notice that the nonlinearity no longer appears in the domain V but directly in 
the Hamiltonian functional. So we derive the linear operator which approximates H^^u) from the 
quadratic form which approximates E{^i^) and obtain the following result. 

Proposition II. 2. The Hessian E"{^i^) of the functional E can he represented as E" {^^){h,k) = 
{Ha, Lin h,k), whcrc Ha, Lin is the operator given by 



H 



a, Lin 



Ha, 

Han 



10 



and 



ai = -iy{2a + 



|2(T 



2(7 + 1 

4tt 



UJ, 



"2 = -V[ 



|2(T 



47r ' 



(17) 



Ha,iin 'is self adjoint with respect to the real scalar product in L (]R'^,C) . 

Proof. The first Gateaux derivative of E{u) reads 
d 



E'{u)[h] = ^{E{u + eh)},=o = ^ [ VMx)-^Mx)dx-i^\qu\^''^iquqh) yu,hGV, (18) 
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while the second Gateaux derivative at reads 

-— {^(^^ + e/i + A/c)K=o,A=o = 3^ / V(l)h{x)-VM^dx 
Tlie last term gives, after some calculation, the contribution (here h = (/ii,/i2), k = {ki,k2)) 



0,A=0 



-{- 



|2(T+2 



e=0,A=0 



dedX' 2(7 + 2'^"' 

So £"'($tj) is given by the direct sum of two quadratic form: one is acting on the real part of the 
functions h and k, and the other on the imaginary part. The term related to the real part is a lower 
bounded quadratic form whose corresponding selfadjoint operator is Ha^^, while the quadratic form 
related to the imaginary part corresponds to the operator (ai and 02 have been defined in pT|) ). 
Then, the operator Ha^un that represents the entire quadratic form £'"($^) is self-adjoint and the 
proof is complete. □ 

Now, to get the linearized equation set u{t) = e*'^*($^ + R) and obtain 

j^R = J{E'{^^) + uj^^) + J{e"{^^) + uj)R + higher order terms ~ J{Ha,Lin + uj)R . 

Summing up, the linearized equation ([6]) becomes 

dR 



dt 



JDR, 



(19) 



where D 



Li 
L2 



, with Lj = + oj, j = 1, 2. Notice that the operator 



JD := L 



L2 
-Li 



(20) 
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is not selfadjoint nor skew adjoint. Nevertheless, a standard application of Hille-Yosida theorem and 
a simple analysis of the resolvent of L which takes into account the factorized structure L = JD with 
D s.a. shows that it generates a semigroup of linear operators with (at most) exponential growth in 
time. A more precise analysis of the resolvent of operator L will be given in Theorem IIV.2I and in 
the appendix O we will prove that the semigroup has in fact a linear growth (see Theorem IC.lh in 
the case here interesting, i.e. a G (0, l/\/2). 



III. ORBITAL STABILITY 



In order to prove the orbital stability of the stationary solutions to equation ([6]), we apply Grillakis- 



20]. As a first step, we recall the following 



Shatah-Strauss theory, and in particular Theorem 2 in 
known fact proved in [3|. 

Proposition III.l. If a{u) = a where a < is a constant, then 

ct(/7„) = {-(47ra)2}u [0,+oo). (21) 

Thanks to the last proposition one can prove the following lemma which implies spectral properties 
needed to verify Assumption 3 in [2^. 

Lemma III. 2. The spectrum of the operator D is 

a{D) = {-4(t((t + l)uj, 0} U [u, +00), 



and ker(Z)) = span 







Proof. Since D is the direct sum of the operators Li and L2 acting on L^(M'^) ©L^(M'^), its spectrum 
is given by the union of cr{Li) and a{L2). From (I2ip follows 

a(i7„J = {-(2fT + l)2w}U[0,+oo), a{Ha^) = {-w} U [0, +00). 

Then 

o-(Li) = o-(FaJ + w = {-Aa{a + l)a;} U [w, +00), (j{L2) = (t{Ho,^) + o; = {0} U [w, +00). 
Hence, ker(Li) = {0} and ker(L2) = spanj^^j}, which concludes the proof. □ 
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We can now prove the following 



Theorem III. 3. (Orbital stability) For each w > 0, i/0 < fi < 1, then <I>^ is orbitally stable. If 
a > 1, then is orbitally unstable. 



Proof. Well posedness and existence of a branch of standing waves, i.e. Assumptions 1 and 2 in 
are proved in fl and 0] and in the previous section, while Assumption 3 is true thanks to Lemma 
IIII.2[ Hence, from Theorem 3 in [2^ we have orbital stability if d"{uj) > and orbital instability if 
d"{uj) < 0, where d"{u}) = ■£;;\\^ujix)\\\2- In order to inspect the sign of d"{uj), we compute 

\^Ux)\\i2 



hence d"{uj) = §^(4^^,)!/^ ^2^*-^^°") which concludes the proof. □ 



A. The case (7 = 1 



Since theorem 3 in [20] does not give information about orbital stability of the stationary state 
gia;t(j)^ when d"{u)) = 0, we need to inspect the case a = 1 apart. In such case, equation ([6]) exhibits 
one additional symmetry (see Q])- 

Remark III. 4. Equation ([6]) is invariant under the pseudoconformal transformation 

In [l| it is proved that equation ([6]) may have some non global solutions which blow up, in the 
following sense: the solution u{t) of equation Q blows up (in the future) at time T < +00 if 

limsup ||V(/)m||l2 = +00. 

Due to the conservation of the energy this condition is equivalent to lim sup^.^.-p- |q'u(t)| = +00. 
Thanks to the pseudoconformal invariance we prove that in any neighbourhood (in energy norm) of 
each standing wave there are initial data of a blow up solution. 

Theorem III. 5. Fix a = 1 and to > 0. For any 6 > there exists a blow up solution u{t) G V such 
that ||n(0) - ^uj\\v < 5- 
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Proof. Applying the pseudoconformal transformation to the sohtary wave e*"*$^ one gets that for 
any T > 0, the function 



_s_ CO 



1 /A _ ^1^1 , 

1/4 g T-t ■ 

e MT-t) 



yJ^KV 47ry/T - t 



X 



'1/4 



e T 



is a solution to equation Thus, for any T > 0, the initial datum ut{x) = e* 

\/A-Kv ATr\/T — t\x\ 

gives rise to a solution that blows up at time T. Now, let uj depend on T as w = wT^, so that 

ut{x) = e * 4t 

We prove the theorem by showing that ||(e~* 4t — — t- as T — t- +oo. Indeed, noting that the 

, |,|2 

function (e"**^ — 1)$^ belongs to 



e-^4T - l)^^\\y = ||V((e-*4T - 1)^^)\\l2 < — III • l^u^h^ + 77^111 • l'V$^||i2 ^0, T ^ +00. 



2T 



4r' 



□ 



IV. SPECTRAL AND DISPERSIVE PROPERTIES OF LINEARIZATION L 



In Section III CI we proved that the linearization of ([6]) around e^^^^ujix) has the structure 



^ = JDR = LR. 
dt 



where we denote 



L = JD 



L2 
-Li 



(22) 



The generalized kernel of L is defined as Ng{L) = ker(L'^). In what follows let us denote 

ken 



/4^' ' 47r ' 



1 

WT: 



4(Taj y47r;^y 47r|x| \^ 47n/ y 47r ' 



~r o,., ,1^ ~r n /-rl'i' 



In Appendix |A| we prove the following theorem. 



47r + 3 1^1 47r 



Theorem IV. 1. If a ^ 1, then Ng{L) = span 

f / 

Moreover, if a = 1, then Ng{L) = span 





M 








\^ 1 









5c. 
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In the following section we provide an explicit description of the spectrum of the non-selfadjoint 
operator L and the dispersive estimates for the action of the propagator e~^* upon the absolutely 
continuous subspace. 



A. The resolvent and the spectrum of the linearized operator 



The purpose of this section is to prove an explicit formula for the resolvent of the linearized 
operator L. For later convenience, it is useful to denote 



Gui±i\{x) 

with the prescription 9-v/^a7±T\ > 0. 



47r|x| 



w > 0, A G 



(23) 



Theorem IV. 2. The resolvent R{\) = {L — XI) ^ of the operator L defined in (I22p is given by 



R{X) 



—XGx2* 


-r^2* 


An 


Ai 


iT.2 


r^2* 


—XGx'2* 


+ Z)(A2)' 


— iSi 


A2 



(24) 



for a// A G C \ {A G C : -^(A^) = 0, or 5i(A) = and |9(A)| > uj}, where the symbol * denotes the 
convolution and 

Gx^ix) = TTT (Gcu-ixix) - G^+ix{x)) , r;^2(x) = ^ {G^_ix{x) + G^+ixix)) . 
2iX 2 

Moreover, 

D{X^) = 327r2aia2 — 4i7r(ai + Q2) ^V— w + iX + V— w — iX^ — 2\/ —uj + iXV—uj — iX, 
and the entries of the second matrix are finite rank operators given by 

Aif = [iA(47ra2 - iV-uj + iX){Qx2, f)L'2 - (47rai - iV-u) + iX)(Tx2, f)L2]Gui+ix+ 



+ [fA(47ra2 - iV-uj - iX){Qx2j)L2 + (4vrai - i\/ -u - iX){T ^2, f)LAGuj-ix, 
Mf = [iX{iiTai - iV-uj + iX){gx2,f)L2 - (47ra2 - iV-uj + iX){Tx2, f)L2]G^+ix+ 
+ [iA(47rai - - iX){Qx2,f)L2 + (47ra2 - iV-uj - iX){Tx2, f)L2]G^-ix, 
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Si/ = -[?A(47ra2 - iV-uj + iXjiGx^, f)L'2 - (47rai - iV-uj + iX){Tx2 , f)L2]Gi^+ix+ 
+ [iA(47ra2 - iV-w - i\){G\2,f )L2 + (47rai - iy/ -uj - iA)(r;^2, /)i2]Gt^_iA, 

S2/ = -[iA(47rai - zV-w + iA)(^;,2, 7)^2 - (47ra2 - i\/ -u + i\){T x^, f)LAGw+iX+ 
+[zA(47rai - W -uj - i\){Qx'iJ)L'^ + (4^2 - W -u) - i\){Vxi, f)LAGuj-i\, 



for all f eL'^i 
Furthermore, the spectrum of the operator L is given by 

a{L) = aess{L) U ad{L), (25) 

where the essential spectrum is 

(Jess{L) = C+ U C_ = {A G C : 5R(A) = and'^iX) > w} U {A G C : 5R(A) = and'^{\) < -uj}, 
and the discrete spectrum depends on the paremeter a as follows: 

(a) if a G (0, l/^/2], then the only eigenvalue of L is with algebraic multiplicity 2. 

(b) if a ^ (l/\/2, 1), then L has two simple eigenvalues :ki2uy/l — u'^oj and the eigenvalue with 
algebraic multiplicity 2. 

(c) if a = 1, then the only eigenvalue of L is with algebraic multiplicity 4. 



(d) if a E (l,+oo), then L has two simple eigenvalues ilas/a'^ — Ico and the eigenvalue with 
algebraic multiplicity 2. 

Proof. Step 1: The Green's function of the operator "H = + (— A + o;)^ 
As a first step, let us compute the Green's function of the fourth order operator 

Lemma IV.3. For any fi E C, oj > 0, the Green's function of the operator % is 

= (G^_i^(x) - , (26) 
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Proof. By definition of Green's function, solves the equation + (— A + tj)^]t/^(x) = 5{x). Taking 
the Fourier transform, one gets 

where the function G^±i^ was defined in ([23|) . The proof is complete. □ 
Remark IV.4. G D{-A) and qg^ = 0. 
Let us define 



■■= (-A + a;)g^(x) = = - \^G^_i^{x) + G^+i^{x] 

Step 2: The resolvent of the operators L2L1 and L1L2 

Exploiting the result of the last step, we can compute the resolvent of the operator 

L2L1 = {Ha2 + uj){Ha^ + w), 

whose domain is 

D{L2Li) = {u G L2(m3) : u £ D{H^^), {H^^ + uj)u G D{H^^)] = 

u G L2(m3) . ^ = ^ + + 5t^' with ^ G Ae G ^2(^3)^ 

ovr 47r|x| 

p,gGC, e(0) = aig, (-A + L^)C(0) =02(^6^ -p)}. 

Remark IV. 5. For later convenience, let us note that the domain of the operator .^2-^1 can be 
written in the following form 

D{L2Li) = {tx G l2(m3) : u = i^+ pG^+,^ + with G ^^(M^), p, g G C, 

= ai{p + q), 



, s ^/— a; — i.fn, a; + i.fa | 

(- A + uj)i^{Q) + ^p:^^—-^ - = a2i^{q - p)\ . 

This fact can be seen by adding and subtracting a suitable linear combination of Gj^+j^ and G^^i^ 
to a generic function of the domain D{L2Li), namely 
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and inspecting the relationship between the boundary conditions one gets the claim. 
Let us denote 

^21 = M + L2L1. 
In the following lemma the inverse operator of 7i2i is constructed. 
Lemma IV.6. For each fi EC, let us define the operator ■ L^(M^) — >■ D{y,2i) as 

where p,q : L^(M^) ->■ C are 



P = i^Din) [^V/^(47ra2 - ^1/^a^+T7M)(^M'/)L2 - (47rQ;i - i^-o; + iy7Z)(r^, 7)^2], 
9 = i^Di^i) [^VM(47ra2 - i^/-<^ - i^){Qt,, f)L^ + (47rQ!i - z^-o; - z^)(r^, 7)^2], 

= 2(47r)^aia2 — 4z7r(Q:i + 02) ^y^-oT+i^ + ^J^-uj^-U/^ - + —uj — 

Then 

^21^21' = 

Proof. First of all let us observe 

Prom the previous equality and the fact that G ff^(M^) follows 

^21^2"// = (/^ + (^0 + * / + pGo^+ivn + QGu^-i^) = 

= + (_A + u;f){f, + (-A + ufr'f = / 
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for all / G L2( 

It is left to determine p and q such that "HgiV ^ D{%2i) = D{L2Li) for all / G L'^{M?). From the 
boundary conditions in the definition of the domain it follows that p and q must satisfy 



{Q,, fh^ + ^P^^^ + ^^^^^ = + ^z) 

^ (r^,/)L2 + '^-^vg _ ^ a2i^{q-p). 

The last system is solved by the functionals p and given in the statement of the lemma. 

To conclude the proof one has to show T-L2i'H2iu = u for any u £ D{'H2i)- For this purpose let us 

set u = + ciG^j^i^ + bG^^i^ for some G H^(U.^) and a, 6 G C, then 

'H^lU2iu = ^2iH/^ + (^0 + + aG^+*v^ + bG^^,^) = 



= 1^21 {i^ + (-A + ujf)if, = Cti + aG^+i^ + hG^_i^ = u, 
since the functionals p and q are defined in such a way that u G D{L2Li). □ 
Remark IV. 7. The resolvent of L1L2 is obtained by switching ai into 02 and vice versa. 



Step 3: Concluding the proof of the resolvent formula ([2? 
Finally, we have all the tools needed to prove formula (j24p . By direct calculations one can show 



R(X) = {L- 



-A(A2 + L2Li)-i -L2(A2 + LiL2)-i 
Li(A2 + L2Li)-i -A(A2 + LiL2)~i 



with A in some open subset of C. This identity is proved in Appendix [Bl 

In order to conclude the proof one can replace A = ^/JI and the expression for p and q in (j27p and find 
the operators Ai and A2. Thus, it suffices to observe that the operators Si and S2 can be obtained 
applying Li and L2 to Ai and A2 respectively. 
Step 4: The spectrum of the operator L 

The statement about the essential spectrum of L is a consequence of Weyl's theorem (see for example 
Theorem XIII. 4 in [3]). On the other hand, the eigenvalues of L are given by the poles of the 
resolvent ()24|) . or equivalently by the complex roots of the function D{X); these can be computed 
through a lengthy but elementary calculation, here omitted. The same analysis gives the presence 
of a resonance at the endpoints of essential spectrum for the case = ^ 
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B. Dispersive estimates for the linearized problem in the case a G (0, l/\/2) 



In this section we focus on the behaviour for large t of the propagator e~^* restricted to the 
essential spectrum and for a G (0, l/\/2)- 

In order to achieve an effective estimate we need to set up a suitable functional space. In our case 
the following weighted spaces are needed 

L^(M3) = |/ : m3 ^ C : j ^ w{x)\f{x)\dx < +oo| , 

and 

L~_i(M3) = // : m3 ^ C : esssup{w{x))-^\f{x)\ < +oo| , 

where w{x) = The use of such weighted spaces seems to be unavoidable, due to the singularities 

of the elements in the domain of the delta interaction (see [ld| for dispersion in the case of delta 
interactions in M^) and the related, more complicated, operators constructed from it and here studied. 



Theorem IV. 8. There exists a constant C > such that 

^ [ [ {R{X + 0) - R{X - 0))e-^*/(y) dXdy < C f 1 + ^"j t'l [ + \fiy)\dy, 

for any f € ^^(M^), where 

C+ = {A G C : K(A) = and 9(A) >uj}, C_ = {A G C : ^{X) = and 9(A) < -u} . 

Proof. One can compute the propagator e~^* as the inverse Laplace transform of the resolvent of L. 
In particular, by Theorem IIV.2I and applying the residue theorem, it follows that for t > 



— [ R{X)e-^^dX = — [ R{X)e-^^dX + — [ {R{X + 0) - R{X - 0))e-^^dX, 
with r G (0,a;) and R{X ± 0) = lim,^o+ R{X ± e)- 

We show the computations only for the component i?i^i(A) of the resolvent whose analytic expression 
is given in (|24p . since the other components can be handled in the same way. 

Recalling the definition of ai and 02 given in equation (IITD . i?i,i(A) can be written as an integral 
kernel as follows 



-1 i\/ —uj+i\\x—y\ i\/ —u)—i\\x—y\ 

RiAX;x,y) = + (28) 
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87r|x||y|[(2(T + l)a; + i((T + l)v^ (-v/— a; — iX + \/—uj + zA) — V— w — zA\/— + 



[{a + 1) + iV-^^ + «A]e^v'-'^+^^(l^l+lf I) + ^^eW-u^+iMyle^V-'^-iXM 
87r|x||y|[(2(j + l)u + i{a + l)-v/w (V— w — i\ + \/—uj + iA) — \/— a; — iX^J—uj + iA] 

Since from equation (j28p it is clear that the computation of the integral on C+ and on C_ are 
exactly the same, we concentrate on the cut C+ only. On C+, —uj + iX is continuous while, by the 



prescription '^{y/—uj it iX) > 0, 

lim \/ —uj — iiX + e) = — lim \/ —uj — i(X — e) = —\/—uj — iX. 
Let us perform the change of variables k = \/ —lo — iX, then 

/• r+oa 

/ (i?i,i(A + 0) - i?i,i(A - 0))e'^^dX = ie-*^* / F{k)2ke-'^^ dk, 

Jc+ ' ' J~oo 

where F{k) is the function i?(A + 0) — R{X — 0) in the new variable. 

The function Ri^i defined in ()28p is the sum of a convolution summand and a multiplication 

summand Rm 1 1 , where 



-uj+iX\x-y\ _ ^i^/-uj-iX\x-y\ 

R*,i,ii>^;x,y) = — ■ 

^ 87r|x — y\ 

and 

aV^e*^-"-^^ls'le^^-'^+*^l^l + \(a + 1) - W-u + iAle^^-^-'^d^l+I^D 

Rm,i,i{X;x,y) = i — + 

87r|x||y|[(2(j + 1)lo + i{a + l)v^ (v — w — iA + y— a; + iA) — y— a; — iX\/—uj + iX] 

[(cr + 1)^ + iV-w + iX\e'^-^+iHM+\y\) + aV^e^^-'^+^^ls'le*^-'^-*^!^! 
87r|x||y|[(2(j + 1)lo + i((T + l)-v/tJ (V— w — iA + \/— a; + iA) — \/— a; — iA\/— a; + iX] 

Then we can define 

F,{k) = i?,,i,i(A + 0) - i?*,i,i(A - 0), and F^(fe) = i?^,i,i(A + 0) - - 0). 

sm{\x~y\k) rp^,,^ Wrvnulo Q Rr^-] -.r. 



One can easily compute F^, and get F^:{k) = . Thus, by formula 3.851 in 

F,{k)2ke-''^'dk = -llif-le*^, 
for any f > 0. Hence, there exists some C > such that 

^. [ F,{k;y)dkf{y)dy < -^t"^ / \f{y)\dy. (29) 
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Let us estimate Fjn{k)2ke~'^^^^ dk. One can notice that Fjn{k) is the sum of terms of the form 
87r | x || j/ |^=*=^(^)^^'^^' where g±s{k) is a rational function of -^—00 — iX and sZ—co + iX multiphed by 
giV-2w-A;^s^ and s can be 0, \x\, \y\ or |x| + Let us consider the term 



g-(\x\+\y\){k)e 



-ik{\x\+\y\) 



{a + l)^/oJ + iV-2uj - A;2 



-ifc(|a;|+|y|) 



(2(7 + 1)0; + i{a + l)^{V-2u; - P - A;) + A;V-2a; - P 

where G C^([0,+oo)) and \g-(\x\+\y\)ik)\ ~ 1^ as A; ^ +00, 

hence G L2((— 00, +00)). Moreover, 

9'-{\x\ + \y\)ik) = 



(2(7 + 1)0; + i{a + 1) ^^^(^-20; - P - A;) + A;V-2a; - A;^ 



(cr + 1)^ + fV-2a; - A;2 



[(2(7 + l)a; + i{a + 1)^/uj{V-2oj - k"^ - k) + k^/-2u - k'^Y 



^ i((7 + 1) + \/-2a; - A;2 



which belongs to L2((— 00, +00)). Since is an element of i7-'-((— 00, +00)), G 

L^((— 00, +00)), where is the inverse Fourier transform of Furthermore, one can 

compute the inverse Fourier transform of 2^6"***^^ in (—00, +00): 



2ke-'^^\-'^'dk 



1 -A 
-e 4it . 



{Ami) 



From the last equality it follows 

r<+00 



-/ / 

27rz 7k3 7_, 



87r|x||y| 



5-(|.|+|2,|)(A:)e-^'=(l-l+l^l)2A:e-**'='(iA:/(y)(i2/ 



1 



+00 

K3 7_oo 87r|2:||y| 



^-(|x|+|2/|)(«)^^t(^i- \A - \y\)duf{y)dy 



7m3 |y| 



dy, (30) 



where the last inequality follows from Holder inequality and C > 0. The other terms in Fm{k) are 
handled in the same way. Summing up, let / G Li,(M3) then 



/ / {R{X + 0) - i?(A - 0))e->^'dXf{y)dy 



< 
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< 



27r 



{R{\ + 0) - R{X - 0))e"^*dA/(y) 



dy+ 



+ f I {R{\ + 0) - R{\ - 0))e-^'dXfiy) dy) =^{I + II). 

JR^ JC- J 

Let us estimate the integral /. Thanks to the estimates (j29p and ()30p one has 

f +00 



F{k)2ke~'"' dkf{y) 



dy< 



< 



+00 



F^{k)2ke-'*^ dkf{y) 



dy + 



Fmik)2ke-'"' dkf{y) 



dy< 



\f{y)\dy + 



\A Jr3 \y\ ) V F 



7r3 V \y\ 



dy. 



The integral // can be estimated in the same way, which completes the proof. □ 

Remark IV.9. Evaluating the propagator e~^* at t = one gets 

1 = — [ R{X)dX + — I {R{X + 0) - R{X - 0))dX = Pq + Pc- 
27rz J\x\=r 27r? Jc+uC- 

From Lemma IV.ll it will follow that the operators Pq and Pc are symplectic projectors onto the 
subspaces associated to generalized kernel and to the continuous spectrum respectively. Finally, let 
us note that explicitely integrating the resolvent around its poles it turns out that the dynamics 
along the generalized kernel grows linearly in time. This fact is proved in Appendix [Cl 



V. MODULATION EQUATIONS 



In this section we restrict to the case cr G (0, l/\/2), summarize the main technical steps and give 
some preliminary results towards the proof of asymptotic stability of standing waves. In particular, 
we write the so-called modulation equations that rule the evolution of a perturbed standing wave 
when splitted in a solitary component and a fluctuating one. We recall once more that the scalar 
product we adopt is the real scalar product on the Hilbert space (L^(R^, C) defined in (jlOp . In order 
to make the reading easier, let us give a, 
the roadmap of the classical papers 



jrief outline of the strategy to be employed. We follow 



also adopted for the model with concentrated 



nonlinearity in dimension one in [4| and [22|. More specifically, we decompose the dynamics in 
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the neighbourhood of the solitary manifold in a "longitudinal" and a "transversal" component with 
respect to the generalized kernel Ng(L), given in Theorem lIV.il of the linearized operator L. In order 
to perform the required analysis, we exploit the symplectic structure introduced in Section III A[ Let 
us begin by noticing that the solitary manifold Ai defined in (I16p is a symplectic submanifold of 
(L^(M'^, C), il), invariant under the flow of ([6]). Its tangent space at the standing wave is two- 
dimensional and is generated by the vectors ■^{e'^^^i^}g=o and ■^{^uj}uj=o, in real representation 
given by 




— {e*^$^}e=o ^ ei = and —{^u}^=o ^ €2 = 

de j dw 

where (p^^ = ^^o; was defined in Section [TVl However, when no confusion arises we use the shorthand 
expressions <I>^ and with the meaning of the corresponding real representative vectors (second 
component vanishing). As already remarked the couple of vectors {61,62} is a basis for Ng{L). It 
immediately seen that r2(6i, 62) = ^^II^^IP 7^ 0, thanks to the condition a G (0, l/\/2) guaranteeing 
orbital stability. So the symplectic form is nondegenerate on the solitary manifold which is a 
symplectic submanifold. By its very definition, M is invariant for the flow of ([6]). 
The following lemma establishes the relation between the spectral projection Pq introduced in Remark 
IIV.9I and the symplectic projection onto the solitary manifold. 

Lemma V.l. Let A = |^||<^aj|li2, then for any f G L^(M^) 

Pof = if, ^u.) J^u. - (/, J^uj) , (31) 



where f2(-,-) was defined in p^ . 

Proof. The explicit expression of the spectral projection Pq = i\\\=r R{^)dX can be recovered by 
Appendix [Cl and the equivalence with the r.h.s. follows by straightfoward calculations. □ 

Notice that the given representation of Pq is well defined thanks to the fact that A > 0, again as 
a consequence of the choice a G (0, l/\/2)- Moreover, Pq is a symplectically orthogonal projection, 
in the sense that given a couple {C)/} with C, G Im Pq and / G Ker Pq) one has f) = . 
In particular it is useful to note that due to the definition of symplectic form 0, a state / with 
vanishing component along the continuous spectrum of L is orthogonal to the vectors Jei and J62, 
or in complex notation, to <I>^ and i-^j^iu = i^uj- 
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After these preliminaries, as anticipated in formula ([S]), we write the solution to ^ as 

^x(t,x) = e*®W(<f^t)(x) + x(t,x)), Q{t)= fu;{s)ds + j{t), (32) 



with the final goal of proving that the solution decomposes in the sum of a solitary component and 
a dispersive one. 

The local splitting of the invariant symplectic manifold (L^(]R^, C), il) in two symplectically or- 
thogonal manifolds, the finite dimensional solitary manifold Ai and the infinite dimensional range 
of the spectral projection on the continuous spectrum suggests to symplectically project the flow 
according to this decomposition (see also Remark lIV.Op . in order to obtain the so called modulation 
equations. The projection along Ai ("longitudinal") gives rise to two ordinary differential equations 
for the frequency to and the phase 7 of the solitary wave, depending parametrically on the fluctuating 
component x; while the projection on the continuous spectrum ("transversal") gives a partial differ- 
ential equation for the remainder x (with coefficients depending on 7 and w). The solution to the 
equation for the x component will be shown to decay in time in suitable norms. As a consequence, 
one has the asymptotic behaviour of the solutions for the parameters to and 7 of the solitary wave, 
to be shown in Section 6, and finally asymptotic stability, which will be the subject of Section 7. 
To deduce the modulation equations it proves convenient to make use of the variational formulation 
of equation ([6]) 

(i^it),v^^^=E'[uimv) yvGV. (33) 



To begin with, we replace in the previous equation the Ansatz (I32p . 

By equation and Proposition 111.21 equation ([33]) can be rephrased as 

(i^(t),t') ^= Qa,Linixit),v) + W)i^At) + Xit),v)L2 + dj{t) (^-z'^^^^v^ ^ + N {q^{t) , {t)) 

(34) 

for any v £ V . 

Here Qa Lin is the quadratic form of the operator D defined in (jl9p and acting as 



Qa,Lin{x, v) = (V(/>x, V0^)i2 - -^'Siiq^q^) - a^^R.q^Rqy + uj{x, v)l2 , 
and the nonlinear remainder N(q^,qy) is given by 

Niqx, Qv) = -z^kx + Q^I^'^^HQx + + ^(2f^ + '^)\quj\'^''^qx^Qv + i^lqu^l'^'^^q^'^qv + iy\quj\'^''^iqc.q^)- 
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Remark V.2. The remainder N{q^,q^) depends nonlinearly on x (and a;) and it is real linear in v; 



such that N{q^, q^) = ^N{q^)'q:^. It is a peculiarity of this model that in fact it depends just on the 
charges of x and v. Moreover, by its very definition the remainder it is the difference between the 
action of the complete vector field and its linear part at the solitary wave, and so it is quadratic in 
Qx near x = 0. 

Corresponding expressions can be given with obvious modification in purely real form, which we 
omit for the sake of brevity. Since uj, 7 and x are all unknown the Ansatz (j32p makes the problem 
under deter mined, and a supplementary condition is needed to give a unique representation of the 
solution; a way to close the system for oj, 7 and x is to require that the x component is decoupled from 
the discrete spectrum, i.e. PqX = 0) or equivalently to project equation (j3l|) onto the symplectically 
orthogonal complement of the generalized kernel of L. The corresponding modulation equations take 
different forms according to the way one writes the projection and we give two of them for future 
reference. In the following we denote hy Ql the bilinear form associated to the linear nonself adjoint 
operator L. 

Theorem V.3. (Modulation equations I) Let x a solution to equation (|34p such that PQx{t) = 
for all t > 0, and let the functions uj and 7 belong to C^{R); then uj and 7 solve the equations 



so, by Riesz representation theorem and with a slight abuse of notation, there exist a vector N{q^) 



(JiV(gx)9Po*(*-+x)) 



(35) 



UJ 




and 
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(36) 



Proof. We adapt the reasoning in [7] 



Equation is equivalent to 



^($^ +X),v] = QUX, v)) + 7 (Ji'^co + x),v)l2 + ?fi{JNiq^)q;;) 



(37) 



/ L2 



Set V = Pq{^ui + x)'i notice that differentiating in time Pqx = 0, one has 




where expressions such as -^X are computed from the representation given in (|3ip . 
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Moreover, one immediately has the identities 



Ql{x, Pq{^u, + x) = QMx, + x)) = 

and, using PqJ = JPq , 
So one remains with 

(poj^i^u. + x), + x) = '^(J^i^xhp.n^^+x) 



from which the equation for to follows. 



Now let us consider the test function JPo^(^aj + x); ^^id notice the following facts, in which use is 



made of JPq = PqJ. 

j^{<i>co + x),JPQ{^u. + x)) ^ = f^(«'^ + x),J^o('J>- + x)) ^ = (Poj^i'^u + x),JPo{'^. + x)] ^=0; 

Ql{x,JPoj^{^. + x)) = ^ ■ 

It follows from the weak equation (|37p 

7 + X, ^o|('J>a. + X)) = 5R(JiV('Zx)9jPo^(ci.„+x) 

and hence, after substituting the expression of Pq-^{(^^ + x) determined above and cancelation of to 
the equation for 7 follows. This ends the proof. □ 

Two properties of the modulation equations which will be useful in the subsequent analysis are the 
following. 

Corollary V.4. Under the hypotheses of Theorem I V. ^ and if it is known that ||x||li, ^-5 sufficiently 
small, the right hand sides of (j35p and (j36p are smooth and there exists a continuous function 
TZ = TZ{uj, llxllLi,) such that, for any t > 0, 

\djit)\<n\q^it)\^ and \j{t)\ < n\q^{t)\^ 

The proof of the previous result is a consequence of two facts. In the first place (ifuj , u) = 
^^||<I>^|p > by condition a S (0, l/\/2) which gives orbital stability; secondarily, the nonlinear 
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part in (|34p actually depends only on the charges and q^, provided that < c, there exists a 
positive constant C > such that the denominators in ()35p and (|36p are strictly away from zero and 



The second property concerns the compatibility of the orthogonality condition of the fluctuating 
part X with arbitrary choices of initial data. The following lemma assures in fact that orthogonality 
condition Pqx = can be satisfied at the initial time in the neighbourhood of the solitary manifold 
without loss of generality. 

Lemma V.5. Let u{t) G C{W^,V) he a solution to equation ([6]) with u(0) = uq ^ V r\ L]^ and 

assume 

for some > and Oq G M. 

Then, there exists a stationary wave e*^''*!'^^, and Xoi^) with Po('^o)xo = such that uo{x) = 
e*^o + xo{x)) , and ||xo||vnLi = 0{d) as d 0. 



The result is commonly stated as a preliminary step in the analysis of modulation equations (see 
for example j2l| . Q and fl). The proof is an application of the implicit function theorem making 



use again of the condition ^H^l'i^lP 7^ 0; we omit details and refer to the quoted references. As a 
consequence of the previous lemma, in all proofs in the rest of the paper we can assume PqXo = 
where xo = x(0)- 

An equivalent form of the modulation equations for the soliton parameters oj and 7 can be ob- 
tained exploiting the characterization of the condition Pqx = through the (Hilbert) orthogonality 
{x,^ui)l'2 = = (x,^¥'tj)L2- III some respects they are more transparent and we give them making 
use of the complex writing. 

Theorem V.6. (Modulation equations II) Let x be a solution to equation ()34p such that Pox(^) = 
for all t > 0, and let the functions uj and 7 belong to C'^iM); then u and 7 satisfy the equations 



^ ^ ((X,<^^)l2 + {iPu;,^u^)L^)N{x,i^oo) - ix,i^u)L2Nix,^o.) .ggs 
(V'a;,^a;)i2 - (X,'P^)i2 
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Proof. Differentiating in time the orthogonality conditions {x^^ui)l^ = = (x,iv'tj)L2, it easily 
follows that 

So testing the weak equation for x with and ip and taking into account properties of operators 
Li and L2 and orthogonality conditions again, one obtains the system 

The thesis follows solving for cj and 7. 

□ 

Notice that to this second form of modulation equations apply similar remarks to the ones made 
for the first form. In particular, if a priori estimates on smallness of x £^re known, the modulation 
equations are well defined thanks to the condition ^H'l'a;!!^ > , and the analogous of Lemma ED 
holds true. 

VI. TIME DECAY OF WEAK SOLUTIONS 

The goal of this section is to provide the time decay of the transversal component x of the solution 
u (see (31)) to equation ([6]); the result we achieve shows that x is in fact not only a fluctuation, but 
a decaying disperive remainder and it paves the way to the proof of asymptotic stability of standing 
waves, that is given in the next section. To this end we follow the idea developed in 0,0], Q] for the 



standard NLS and applied in [4] to the case of 1-d concentrated nonlinearities. 
For any T > 0, define preliminarily the so-called majorant 

M{T)= sup [(l + t)3/2||;^(t)||^^ +(i + t)3(|^(t)| + |^(t)|)l. (40) 

0<t<T L J 

We aim at proving that the majorant is uniformly bounded in T by a constant M = 0(d) where d 
is the size of the dispersive component x- The proof of such bound is the content of the following 
theorem. 
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Theorem VI. 1. Let u € C(]R"^, V) be a solution to equation ^ with u{0) = uq £ VnLj^ and define 
d := ||uo - e^^o^^JI 

ynLi ) for some ooq > and 9q € M. Then, if d is sufficiently small, there are 
(^,7 G C^(M^) which satisfy (|35p -(j36 p . and such that the solution u can be written as in ()32p . 
Moreover, there is a positive constant M > 0, depending only on the initial data, such that, for any 
T >0, one has M{T) < M, and M = 0{d) as d-^0. In particular 

||x(i)l|L- , <M (l + t)-3/2 vt>0 (41) 

\j{t)\ + \dj{t)\ <M (l + t)-^ yt>0. (42) 

The previous theorem is imphed by the following proposition that is proven in Section IVI CI by using 
the results given in Sections IVI Al and lVI Bl and the dispersive properties of the linearization operator 
L given in Section [IV Bl 

Proposition VI. 2. Under the hypotheses of the previous theorem, assume that there exist some 
ti > and p > such that M(ti) < p. Then there are two positive numbers di and pi, independent 
ofti, such that if d = ||xo||ynLi, < ^^i ^'^'^ P < Pi; then M{ti) < |. 

Indeed, if Proposition IVI.2] were true, then Theorem IVI. II would follow from the next argument: let 
I C [0, +oo) be defined as 

X = {ti > : 7 G C\[0, ti]), M{ti) < p}. 

I is obviously relatively closed. On the other hand, the thesis of Proposition IVI.21 and the estimates 
of Corollary IV. 41 imply that X is also relatively open. Hence, the uniform estimate of Theorem IVI. II 
follows from the fact that supX = +oo. 

A. Frozen linearized problem 

Note that the equation ()34p is non autonomous. In order to make its study simpler, it is useful to 
exploit a further reparametrization of the solution x(^)- We fix a time ti > and denote ui = uj{ti) 
and 7i = 7(ii). Now define (in vector notation; we recall that J corresponds to — i) 

e"-^®^*^x(t,a;) = e--^®Wr?(t,x) , where e(t) = + 71 . (43) 

The function r] satisfies the equation 
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dt 



L2 



■ UJ 



doj 



L2 



We need a further manipulation which allows to rewrite the previous equation in a form which 
makes the role of reparametrization clear. To this end we need the following identities, which can 
be obtained from straightforward computations 



jgJ(0-0) ^ e^(®-®) J; 

g^(gj(e-e)^^^) - e-^(®-®)QL(n,z;) = sin(e - @)asqu7h, for any u,v€V, where 



27r 



c^3 



1 
-1 



Making use of the previous identities, one rewrites the equation for r] as 



UJ- 



doj 



,f) + (44) 

L2 



+e-^(®-®)^^±i^sin(e - e)(T3g,g;; + e-^(®-®)jiV(e^(®-®)(?,)g;;, 'iveV . 
Let us define the linearization frozen at time ti as L/ = L(a;i), and observe that for all v € y 



Ql{u.,v) - Qlj{u,v) 



47r 



TguQ^ - (wi - uj){Ju,v)l2, 



where T 



-1 
2a + 1 



Hence, equation (j44p becomes 



dfj 



(45) 



L2 



where the time dependent nonlinear remainder (including now "dragging" terms due to reparametriza- 
tion) is given by 



7J4>^ - w-j— ) ,f ) H j^^^^ Tqr^q^ 



doj 



L2 



+ e 



-J(0-0) 



6) (f^ + 1)\/^ 
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sin(G - Q)azqr,qv + e 



47r 



-j(0-e)j^^gj(0-e) 



(46) 
(47) 
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The gain in changing from original (|34p for the dispersive component to equation (|45|) is that the 
latter is still non autonomous, but now the generator of the evolution is (in weak form) a sum of a 
fixed linear vector field (the frozen linearization Lj) and a nonlinear time dependent perturbation 
(see also This allows to use the known dispersive properties of linearization operator L described 



B. Duhamel's representation 

In this subsection we write the equation (I45|) in Duhamel's representation, to better exploit the 
dispersive properties of the propagator e^'''. This is not a completely trivial task since our frozen 
equation is a variational equation and cannot be written in strong form. In order to reach our 
purpose, we consider (j45p separating in the test function v the regular and singular part, accordingly 
to (I5|). So we begin by setting v = (p^ £ H^(M.^). We get 



in llVBl 




where //(t) = e--^(®W-0W) (7(t)J$a;(t) - ^(t) 




) 



. Hence, by Duhamel's principle one gets 




If one considers the same equation with v = q^Gx for any £ C, one has 




where 



gi{t) = e-^(0W-§W) (4^/A(a + l)y^sm{Q{t) - Q(t))asq^{t)Gx+ 



+8^^/AJ7V(e^(®W-®W)(7r,(t))GA) + 2VX{y^ - ^)Tq^it)Gx. 



And hence, 




Summing up, for any v GV, the solution to equation (|45p can be written as 



{rj,v)L2 = (e^^'^vo + f e"^'^'-'^ fi{s)ds,v\ + ( gi{s)ds, q^G, 



A 



) 




In what follows we will use the following estimate on the function gj. 
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Lemma VI. 3. Under the hypotheses of Proposition \VI.S\ there exists a constant C > such that 

\\9i{t)\\vrMl <C(|g,,P + /5|g^|), 

for any t <ti. 

Proof. First of all let us notice that it is possible to chose ti in such a way that a;(t) > c > for any 
<t<ti, then 

\^/^ - ./IJ^l < C\uj{t) - uji\ < C [ ' \uj{s)\ds < C sup [{1 + tf\dj{t)\] [ \l + s)-^ds <Cp, 

Jt 0<t<ti Jt 

and 

\Q{t)-e{t)\< [ r \uj{T)\dTds+ r \^{s)\ds < p f I {l + T)-^dTds + p [ \l + s)-^ds <Cp. 

Jo Js Jt Jo Js Jt 

The result follows since 



□ 

We end the section with a technical result that allows to transfer dispersive estimates on the frozen 
fluctuating component Pc{Li)r] = Pc{u}i)r] into estimates on r]. This is needed because ij appears in 
the integral Duhamel's equation where estimates have to be done, but the dispersive behaviour is at 
our disposal for Pc{uji)rj. This is stated in the following lemma (see for analogous construction, for 
example, Q Q)- 

Lemma VI. 4. Let the hypotheses of Proposition I VI. S\ hold true and suppose that the quantity 

sup {\u;{t) - uji\ + \Qit) - e{t)\) = 6 

0<t<ti 

is sufficiently small; then, for any t G [0,ti] there is a bounded linear operator n(t) : Pc{u;i){V Ci 
L'^_i) — )• V r]L^_i, and a positive constant C = C{6,uji) > such that ri(t) = Il{t)h(t), and 

C{d,ujiy'^\\h\\vnL°°, <\\v\\vnL'=°_, < C{6,uJi)\\h\\vnL°°_, ■ 

Proof. We give only a sketch of the standard proof, referring for details to the literature cited above. 
Set 7]{t,x) = Po{uji)r] + Pc{uJi)r] = iki{t)<^>^^ + k2{t)-£j^^u^^ + h{t,x) . The condition PqX = 
makes time dependent functions ki and k2 to satisfy a linear system with a source term depending 
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on h; the coefficient matrix has an inverse uniformly bounded in t and ti thanks to the conditions 
(^•^^ ^ const > and {^uji, > const > valid for — small enough. This 

gives a representation of ki and k2 in terms of h and as a consequence the required bound on the 
finite dimensional component. Now define Il{t)h{t) = rj{t) — iki^^Ji + k2-^^u)i and the complete 
bound follows. □ 



C. Proof of Proposition rVI.2l 



Estimate of I7I + |w| . 

Lemma VI. 5. If t] £ ynL°^„i then \qJ < 47r||7/||Loo . 

W Til 1 

Proof. Since r] G L'^^i{R^) then ||r?||L-_^ = sup^,gK3 j^\(t>vix) + 



> 



An I 



□ 



From the last lemma and Corollary IV. 41 one gets 

ItWI + I'^WI <ck^(t)|2 <ci||r?(t)||ioo <ci(l + t)-3M(t)2, vte [0,ti], 

with ci independent of ti. Hence, one can choose pi < ^ and get (1 + t)^(|7(t)| + |^(i)|) < cip^ < 
f, VtG[0,ti]. 
Estimate of ||7?||l°o . 

As explained in the previous section, for any t G [0, ii] we have = PQ{u}i)r]{t) + Pc{uji)r]{t) (for 
the definitions of Pq and Pc see Remark IIV.9P and thanks to Lemma IVI.4I we have r](t) = Ilh{t) 
where U{t) : Pc{uJi){V n L'^_^) V n L'^^^ is bounded. 

In order to estimate ||7/||ioo_ we make use of the equation for h. For all v & V, h is a solution to 

hence, for any v £ V, h satisfies 

{h,v)L2 = (^e'^''ho + j\''''^'-'^Pc{u:i)fi{s)ds,v^ ^ + (^j\'''^'''^P,{ujMs)ds,q,G^ 



L2 



In addition let us assume that v & V r\ , hence by Holder inequality 



{h,v)L2 < I \\e^'^ho\\vnL-° , + 



'M'-'^Pc{uji)fi{s)d, 
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+ 



\qvG\\ 



Now we can apply the dispersive estimate proved in Theorem IIV.8I and get 

||e^^*/io||ynL- , < c(l + t)-^^^ho\\vnLi < c(l + ty^^^d. 

Furthermore, 



<c [\l + t-sr'/^fj{s)\\vnLlds< 
Jo 



<cl (l + t-s)-^/\\j{s)\ + \uj{s)\)ds<c I (l + t-s)-^/^||r?(s 



-3/2 1 



ds. 



Finally, 



'M'-'^P,{ujMs)ds 



vnL° 



<c I {l + t-sr'/%j{s)\\yr,L}ds< 



< c [\l + t-s)-^/\Ms)\\loo +p\\rjis)U^ ^)ds. 
Jo 

Let us define 

m{t) = sup (l + s)3/2||r,(s)||i^ 

s&[0,t] 

Now, taking as the norm of the intersection of two sets the maximum of the norm one has 
and hence, exploiting the duality pairing defined by the inner product in L^, it holds 



+ 



3/2 I 



L^_^ < (1 + tf/'^\\v\\v*nL^_, = (1 + tf^'^ sup 



3/2 



L2 



< 



<c(^d + m\t) j\l + tf'^{l + s)-^{l + t - s)-^/^ds + pm{t) j\l + tf'^{l + 5)^=^/^(1 + t - s)-^'^ds 

Observe that the constant c and both integrals appearing in the last inequality are bounded inde- 
pendently of t, and this implies that for any t G [0, ti] we have 

m{t) < c{d + m^{ti) + pm{ti)) < c{d + pj) < C2d, 

provided d and p are small enough. Since the constant C2 does not depend on ti, we can choose 
d < ^ and finally get 

P 



m{ti) < 



concluding the proof of Proposition IVI.2I 
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VII. ASYMPTOTIC STABILITY 



Now we are in the position to prove the asymptotic stabihty result as stated in the next theorem. 
Before formulating the result, let us denote by Ut the integral kernel which defines the propagator 
of the free Laplacian in M^, namely Ut{x) = {Airit) ' e*~. 

Theorem VII. 1. Assume a S (0, l/\/2)- Let u{t) € C(M+,F) he a solution to equation ([6]) with 
'u(O) = uq r\ L\j and denote d = \\uq — e^^^^o^ollynLi,) foi^ some ujq > and 6q G M. Then, if d is 
sufficiently small, the solution u{t) can be decomposed as follows 

u{t) = e'''^'<^^^ + Ut*^oo + roo, (48) 
where ujoo > and ipoo, ^oo £ iF'i^^), with ||roo||L2 = 0(t^^/^) as t ^ +oo. 

Proof. Along the proof we assume that i-b(^o ~ e^^^^^uo) = 0, and we recall from Lemma IV.51 that 
there is no loss of generality in this choice. First of all let us notice that Theorem IVI.ll implies 
oj{t) — bJoo, and G(t) — ujoot — >• 0, as t — +oo. Next, let us define the modulated soliton as 

s{t,x) = e'''^'^^^^t){x), 

and the function 

x) = u{t, x) - s{t, x). (49) 

By equation ([6|) and (fTSj) one has that, for any v £ V, z{t) is also a solution to 

( ^37' ^ ) = K /" Vcpz ■ V(j)vdx - i^K((|gu|^'^gu - {QsI"^" qs)q^) + ('js - iuj^,v) . 
\ at J ^2 Jr3 V duj J ^2 

As one can verify by directly differentiation, the solution of the last equation can be expressed as 

z{t,x) = Ut* zo{x) + i Ut-r{x)q^{T)dT - i Ut-r * f{s{r))dT, (50) 

Jo Jo 

where we denoted f{s) = js — ioj^ and according to (^91) qz{t) = quit) — (?s(t). Let us consider the 
last integral in formula (|5U|) 



Ut-r * /(s(r))dr = Ut* U.r* /(5(r))dr - / Ut-r * f{s{T))dT, 
Jo Jt 

and note that the regularity of s{t,x) implies tpiix) = U-r * f{s{T))dT G L^(]R^), and ri{t,x) = 

— Ut-r * f{s{T))dT G L^(R'^). Moreover, from Theorem IVI.ll and the unitarity of the evolution 

group of the free Laplacian we have ||n(i)||L2 = 0(t~^), t — )• +oo. 
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To conclude the proof it is left to prove a similar asymptotic decomposition for the first integral in 
the formula (1501). As before, one can write 



/ Ut-r{x)q.{r)dT = Ut* U-r{x)qz{T)dT - / 
Jo Jo Jt 



Ut-r{x)qziT)dT. 



First of all one needs to show that 'il^o{x) = U-r{x)qz{T)dT belongs to L^(R^). To this aim, let 
us observe that i^oix) = ^ (t)' '^ith h{y) = eT^^I'^ T~'^l'^qz (r) dr, hence 



,2 1 



(4vr) 



1 



(2vr)2 

From the first and the last terms one gets tpo S L^(M^) if and only if h £ L^(]R"'", y^dy). On the 
other hand, one can perform the change of variables n = ^ in the integral function h and get 

%)= r e-'y''^^qz(-)dn= H e-^y^-qz ( -) V^du, 
Jo \uj Jo u \uj 

where we set y = -^j-. Then h{u) = \qz {\). Moreover, by Theorem |VLT1 I;^^^ (^)| y/u < jf^p^ 
then h G (M"*" , -y/udn) and hence, by Plancherel's identity h G L'^iW^ , yfydy). 
Finally, let us denote tq = Ut-T{x)qz{T)dT. As before, we have ro(x) = g(^^^, with g{y) = 
e~^y^^^~'^\t — T)~^^'^qz (t) dr. Moreover, we can set y = exploit the change of variables 
u = —jhf in order to get 



g{y)= r e-'y^-qz(t+-)v^du. 
Jo u \ uj 



Again, Theorem IVI. II implies that g{u) = ^qz + ^) G L^(M+, ^/udu), for any t > 0. In particular, 

n ,,2 r , . ^-5/2 

for any t > 0, with c, c > independent of time. Summing up, Plancherel's identity allows us to 
conclude ||r"o||j;^2 = 0(t~^/^) as t — ?■ +oo. 

Hence the theorem follows with V'oo = -zo + V'o + V'l ) aiid = ?"o + ^'i • □ 



Appendix A: The generalized kernel of the operator L 

The aim of this appendix is to provide the proof or Theorem llV.li 
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Proof. It is easy to see that c^^^, with c G C, is the unique family of distributional solutions to the 
equation 

—Au + uu = 0. 

Furthermore, belongs to D^H^^) but not to D(Ha^) since the boundary condition is not satisfied. 
Hence 



ker(L) = span 







Let us now consider the operator 



-L2L1 
-L1L2 

Since the operator Li is invertible, the following holds 

u G ker(LiL2) <^ u e ker(L2) <^ ker(LiL2) = {^oj}, 

u G ker(L2-Li) 4^ u e ker(L2) = {^uj} ^ e D{Hai) such thatLitt = $0,. 
Solving the former equation one gets that ker(LiL2) = {tpuj}- From this follows 















span < 


[1 




•1 







The operator has the following form 



—L2L1L2 
L1L2L1 



As before 



u G ker(LiL2-Li) <^ Lm G ker(LiL2) = {^oj} ^ ker(LiL2Li) = {^poj}, 



G ker(L2-LiL2) G ker(L2) = or L2U G ker(L2-Li) = 

Let us notice that the equation 

— Au + LUU = (fi^ 
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has a unique family of distributional solutions given by 



u{x) = 



A-KV 



+ 



ci e 



16(72a;2y 47r|x| 2^ Att\x 



1 e~^l^l 

8a; 47r Vsawi 8a;i 



47r 



Notice that one must impose that u belongs to D{Ha2) which means that u € L'^{W^) and satisfies 
the boundary condition. This is equivalent to ask the following algebraic conditions to be verified 



ci = 



C2 



3"' 



Summing up, if cr 7^ 1, ker(L2LiL2) = {^oj}- Hence 

ker(L3) = ker(L2), 

which concludes the first part of the theorem. 

While, in the case cr = 1 we get ker(L2LiL2) = {^wj^w}) then 







ker(L^) = span < 
With analogous computations one can prove that 

ker(L^) = ker(L^) = span < 
which concludes the proof. 














9uj 



K 




□ 



Appendix B: Proof of the resolvent formula 



In this appendix we prove the formula 



R{\) = {L- xiy^ 



-A(A2 + L2Li)-i -L2(A2 + LiL2 



(-1 



Li(A2 + L2Li)-i -A(A2 + LiL2)-i 
for the resolvent of the linear operator L. More precisely, we prove the following proposition. 

Proposition B.l. // A G C \ a{L), then R{X){L - \I)u = u, Vw G D{L), and {L - XI)R{X)f = f 
for any / G (L2^™>3^^2 
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Before proving the former proposition, let us prove the following lemma. 
Lemma B.2. For any A G C \ (t{L) the following identities hold 

1. (A2 + L2Li)-^L-^ = L-\\^ + LiL2)-\ 

2. (A2 + LiL2)-i = (A2L^i + L2)-iLrS 

3. (A2 + LiL^)-iL^"^ = Z^"\a2 + L^Ll)-^ 

where L2 is the restriction of the operator L2 on the projection of its domain onto the continuous 
spectrum of L2 ■ 

Proof. First of all, let us notice that all the inverse operators are well defined since A is not allowed 
to be a spectral point of L, Li is invertible and L2 is restricted to a subspace on which it is invertible 
too. 

In order to prove [H we prove the following claim 

(A^ + L2Li)-^L:^^ = (A^Li + LiL2Li)-^ = ^^^(A^ + LiL2)~\ 

To this purpose, let us take any ^ £ L^(M^), then one has 

{X^ + L2Li)-^L^^^e D{L2Li) and L'^^ e D{Li). 

Hence, the following chain of equalities holds 

(A^Li + LiL2Li){\^ + L2Li)-^L^^i = Li{\^ + L2Li)L^\x'^ + L2Li)-^L^^^ = LiL^^i = ^. 

On the other hand, let us take rj £ D(LiL2Li), and observe that, in particular, r] € D{L2Li). This 
justifies the following equalities 

(A^ + L2Li)-^L^\X^Li + LiL2Li)ri = 

= (A2 + L2Li)-^L^^Li{\^ + LsLi)?? = (A^ + L2Li)-^{\^ + LaLOr/ = r?, 

which concludes the proof of the first equality of the claim. The second equality is proved in the 
same way. 

The proof of [3j can be done in the same way switching Li with L2 and L2 with Li. 
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It is left to prove H. To do that, let ^ be in L'^{R^), then (A^L^"^ + Ls)"^^^"^^ G /^((A^Lj"^ + L2)) 
and L^^^ G D{Li). Hence, we have 

(A^ + LiL2)(A2l^i + L2)-^LlH = LiiX^L^^ + L2){X^L^^ + L2)-^L^^C = 6 

On the other hand, for any rj G D{LiL2) one has ?7 G L>(L2) C L'^{R^) = D[L^^), which justifies 

(A^L^i + L2)-1L^1(A2 + LiL2)r? = (A^L^^ + L2)-^L^^Li{X^L^^ + L2)?7 = r/. 

□ 

We can now prove the proposition. 

Proof. I step: proof of the first identity. 

Let us recall that for u G D{L) holds 

RiX){L - XI)u = 

-A(A2 + L2Li)-i -L2(A2 + LiL2)-i 
Li(A2 + L2Li)-i -A(A2 + LiL2)-i 

where 

wi = A2(A2 + L2Li)-\i + L2(A2 + LiL2)-^LiUi - A(A2 + L2Li)-^L2U2 + AL2(A2 + LiL2)-in2, 
and 

W2 = X^iX^ + LiL2)-^U2 + Li{X^ + L2Li)-^L2U2 + X{X^ + LiL2y^LiUi - XL^iX^ + L2Li)-\i. 

We will concentrate on the first component wi, because the second one can be treated in the same 
way. 

n 

The spectrum of the selfadjoint operator L2 is ([3|]) 

a{L2) = {0} U [a;,+oo), 

where is a simple eigenvalue and ker(L2) = span{<I>(^}. Hence, any U2 G D[L2) can be decomposed 
as 

U2 = a^^ + 52, 
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-A L2 
-Li -X 





where a G C and g2 belongs to the projection of D[L2) onto the continuous spectrum of L2. 
Moreover, since L2^u) = 0, one gets <I>(^ G D{LiL2) and 

which is equivalent to (A^ + ^1^2)^"'^$^^ G ker(L2). 

As a consequence, since Li and L2 are invertible on their domains, one has 

wi = A2(A2 + L2Li)-^L:^^LiUi + L2{\^ + LiL2)-^LiUi+ 

-X{X^ + Z^Li)"^Z252 + aZ2(A2 + LiL2y^L2^^L2g2, 
hence, by lemma [B.2I it follows 

wi = (A^L^i + L2)(A2 + LiL2)-^LiUi - X{X^ + L2Li)-^L2g2 + AL^L^^^A^ + L^Li)-iZ^^72 

= (A^L^i + L2)(A2l^i + L2)-^L:^^Liui = m. 

Summing up, we proved 

R{X){L - XI)u = u yueD{L). 

II step: proof of the second identity. 

First of all let us recall that for / G (L^(R^))^ one has 

{X'^ + L2Li)-^fi£D{L2Li) and (A^ + L1L2)- V2 € Z)(LiL2). 
Hence, the following identities hold 

(L - A/)i?(A)/ = 



-A L2 
-Li -X 



-A(A2 + L2^i)-^ -L2(A2 + LiL2)-i 
Li(A2 + L2Li)-i -A(A2 + LiL2)-i 



fl 

f2 



(A2 + L2Li)(A2 + L2Li)-i/i 

(A2 + LiL2)(A2 + LiL2)-V2 



/, 



which concludes the proof. 
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Appendix C: The dynamics generated by L along the generalized kernel 



In this appendix we estimate the behaviour of the propagator around the eigenvalue 0. This is 
achieved in the following theorem in which it is proved that the dynamic has a linear growth in time 
along the generalized kernel. 

Theorem C.l. For any r G (0,a;) the following holds 



l-CT ■K\x\\y\ 



1-a 2Tr\x\\y\ 



{2a^\y\ - 1) 



for any x,y G 



Proof. Since the convolution term of the resolvent R{\) is continuous in zero it suffices to compute 
the integral of the multiplication term. First of all, let us note that the function 



/(A) 



-Ai(A)e ^* 



87r|a;||j/| 



-(47rai - iV^u + iA)e'^-^-^^l^l + {Anai - W-u - iA)e'^-^+^^l^l 

87r|x||y| 



87r|x||y| 



2(47ra2 + V^)e-^l^l (^Mg-V^M A + o(A2)^ + STraie"^!^! (^Mg-V^l^l A + o(A2)^ + 



+2ze-^l2'l ( ( _ + |x| ) e-^l^lA + o(A' 



a; 



1 -a 
2a; 



A^ + o(A4) 



— „ I II I [(47rQ!2 + Vw)|y| + (47rai + Vuj)\x\ + 

1 — a 27r\x\\y\ A 

as A —7- 0. Hence the function /(A) has a pole of order one in zero. Then, by the Cauchy theorem 

one gets 



1 



Ai(A)e-^*cZA 



uj e 



-V^(\x\+\y\) 



1 — 0" 27r|x||y| 



-[(47ra2 + v^)|y| + (47rai + A/a;)|x| + 1] 
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^ [-2aV^|x| + 1]. 



I- a 27r|x||y| 
Switching ai to a2 and vice versa, it follows 



2^ i,.. = 2.|.||,| [(^""^ + + ^^^"^ + ^^1^1 + = 



[-2aV^|y| + l]. 



-At 



^ g-V^(|a;| + |2/|) 

"~1-C7 27r|x||y| 
On the other hand, the function 

is the sum of a continuous function and a function with a pole of second order in zero, namely 

<7(A)e-^* = 

_ (47^al-^v/^^J+a)e-^/=^N + (87^al-^V^^:7^)e^^/^ .y^y. 

D(A2)47r|x||y| " 

Note that (7(A) = Ylik=2 ^k^k with 

1 — 0" 7r|x||j/| 



then, by residue theorem. 



27rz L>(A2) ^ 1-C7 7r|x||y| (1 - (7)7r|x||y| 



'|A| 

In the same way 



2m DiX^y 

which concludes the proof. □ 
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